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Abstract 

Missions to Mars will require electric power. A leading candidate for providing power is solar 
power provided by photovoltaic arrays. To design such a power system, detailed information on 
solar radiation availability on the Martian surface is necessary. The variation of the solar radiation 
on the Martian surface is governed by three factors: (1) variation in Mars-Sun distance, (2) variation 
in solar zenith angle due to Martian season and time of day, and (3) dust in the Martian atmosphere. 
A major concern is the dust storms, which fc occur on both local and global scales. However, there 
is still appreciable diffuse sunlight available even at high opacity, so that solar array operation is 
still possible. Typical results for tracking solar collectors are also shown and compared to the fixed 
collectors. During the northern hemisphere spring and summer the insolation is relatively high, 2-5 
kW-hr/m 2 -day, due to the low optical depth of the Martian atmosphere. These seasons, totalling 
a full terrestrial year, are the likely ones during which manned missions will be carried out. 


1 Introduction 

Missions to the Martian surface and activities there will require electric power. Unmanned survey 
missions may require power levels ranging from hundreds of watts to a few kilowatts, a manned 
Mars base would be likely to require tens to hundreds of kilowatts; while a facility to process in- 
situ resources such as rocket propellants could require power on the order of megawatts. A leading 
candidate for providing this power is solar power produced by photovoltaic arrays. To design 
such a power system, detailed information on solar radiation availability on the Martian surface is 
necessary. 

Mars is a challenging environment for the use of solar power. The twelve-hour night, although 
not as severe as the moon’s 14 day darkness, still requires that any solar power system include a 
large storage system. Mars is further from the sun, and thus the solar intensity in Mars orbit is 
less than half that at Earth (the eccentricity of the orbit means that the actual intensity vanes by 
±19% during the year). Furthermore, the atmosphere is filled with dust. Even on clear days the 
sky is partialy obscured with atmospheric dust, and during dust storms the atmosphere becomes 
nearly opaque. 


* Tel-Aviv University, Faculty of Engineering, Ramat-Aviv 69978, Israel. NRC Resident Research Associate at 
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Solar irradiance at the surface of Mars varies with season and with the amount of atmospheric 
dust [1]. A major concern is the dust storms, which occur on both local and global scales, and 
their effect on solar array output. Global storms occur near the perihelion, and may persist for 
long periods of time, such that the requirement for energy storage quickly becomes much too large 
to be practical. However, as we have shown [2, 3], there is still an appreciable diffuse component 
of sunlight available even at high opacity (optical depth), so that some solar array output will still 
be possible. 

From a photovoltaic system design point of view, the intensity, frequency, and duration of these 
storms may be viewed as “partially cloudy” and “cloudy” days for which additional energy storage 
must be taken into account. 

The full sky, or “global,” solar radiation is composed of a direct (unscattered, or “beam”) 
component, and a diffuse (scattered by the atmosphere) component. Given the optical depth of 
the atmosphere [4], the global radiation is calculated from the normalized net solar flux function 
based on calculated multiple scattering of the incident radiation integrated over multiple incident 
wavelengths [5]. The direct beam component is derived from the optical depth of the atmosphere 
using Beer’s law, and the diffuse component is the difference of the global and direct radiation. 


2 Optical Depth 

For our purposes, the intensity of the dust storms on Mars can be defined in terms of the resulting 
opacity of the Martian atmosphere. This is most conveniently measured in terms of the optical 
depth r, defined as the logarithm of the attenuation of the direct component of a beam penetrating 
the atmosphere perpendicular to the surface. An optical depth of 1 corresponds to attenuation of 
the beam by a factor of e. 

The most direct and probably most reliable estimates of opacity are those derived from Viking 
Lander imaging of the Sun. Figs. 1 and 2 show the seasonal variation of the normal incidence of the 
optical depth at the Viking Lander locations VL1 and VL2, respectively. The season is indicated 
by position of Mars in its orbit around the sun, measured by the areocentric longitude of the Sun 
L $y where at the Mars vernal equinox L s — 0°. 

The measurements are at the sites of the two Viking landers. Viking Lander 1 (VL1) is located 
at 22.3° N latitude and 47.9° W longitude, and Viking Lander 2 (VL2) is located at 47.7° N latitude 
and 225.7° W longitude. Figs. 1 and 2, derived from ref. [4] and discretized for each 5°, show the 
optical depth at the two sites. The optical depth is assumed to remain constant throughout the 
day. Opacities are minimum during the northern spring ( L s = 0° to 90°) and summer ( L s — 90° 
to 180°), and maximum during southern spring ( L s — 180° to 270°) and summer ( L $ = 270° to 
360°), the seasons during which most local and major dust storms occur. When dust storms are 
not present, the optical depth is typically about 0.5. Two global dust storms occurred during the 
periods of each observation, as indicated by the high values of the optical depth. These values 
are lower bound values; the actual optical depth may be higher. The two global dust storms that 
were observed by the Viking Landers make 1977 one of the worst years of dust activities on Mars 
observed. It should be noted, however, that there may be some Mars years which have no global 
dust storms [1]. 
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3 Solar Radiation 


3.1 Radiation at the Top of the Mars Atmosphere 


The variation of the solar radiation at the top of the Mars atmosphere is governed by the location 
of Mars in its orbit, and by the solar zenith angle. The solar radiation is entirely direct beam 
radiation. The beam irradiance, in W/m 2 ) is given by [1,2]: 


G ob = 590 


[1 + e cos {L s — 248)] 2 

( 13^)2 


( 1 ) 


where e = 0.093377 is the eccentricity of the orbit of Mars, L s is the areocentric longitude, and 
248° is the areocentric longitude of Mars perihelion. (Here we use the subscript ‘o’ for air mass 
zero, indicating intensity above the atmosphere, and the subscript b to indicate the direct (beam) 
component). This irradiance is shown in Fig. 3. 


3.2 Solar Radiation on the Surface of Mars 

The variation of the solar radiation on the Martian surface is governed by three factors: (1) the 
Mars-Sun distance, which varies the intensity at the top of the atmosphere, as discussed in the 
previous section, (2) the solar zenith angle, z, and (3) the opacity of the Martian atmosphere. The 
global solar irradiance G is the sum of the direct beam irradiance* G b) the diffuse irradiance Gd , 
and the irradiance due to albedo G a \ (i*e. , light reflected to the array from the surface): 

G = G b + G d + G ai (2) 

The direct beam irradiance on the Martian surface normal to the solar rays is related by Beer’s 
law to the optical depth, r, of the intervening atmospheric haze: 

G b = Gob exp [— r/ cos z] (3) 

where z is the sun zenith angle. 


3.2.1 Horizontal surface 

The beam irradiance on a horizontal surface is given by: 


Gbh = G ob cos z exp (~r/ cos z) (4) 

For a horizontal surface G a l = 0, and the global irradiance Gh is: 


G bi — Gob COS Z 


/(*, r, al) 
(1 - al) 


( 5 ) 


where al is the albedo. /(z,r, al) is called the normalized net solar flux function and accounts 
for the light which penetrates the atmospheric dust. The net solar flux integrated over the solar 


3 


spectrum on the Martian surface was calculated by Pollack [5] based on multiple wavelength and 
multiple scattering of the solar radiation. Derived data from this calculation can be approximated 
by a polynomial expression: 


f{z, t , al) = 


±±±«W ■(■&)’ 

t=0 j=0 k=Q X UU/ 


(a/)* 


(1 - al) 


( 6 ) 


where p(ij,k) are the coefficients of the polynomial given in Table I. The diffuse irradiance is 
obtained by subtracting the beam irradiance from the global irradiance (eqs. (4) and (5)). 

The daily global insolation (integrated solar energy per unit area) on a horizontal surface is 
obtained by integrating the irradiance over the period from sunrise to sunset: 


H h = 


24 . 65 * 

27 r 



cos z 


/(z.t-.q 0 

(1 - al) 


du) 


( 7 ) 


where u> sr and u> ss are the sunrise and sunset hour angles, respectively. The daily beam insolation 
on a horizontal surface is obtained by: 


24.65* T“ 

Hbh. - —7 / cos z ex 

ZTT Jw 9 r 

and the daily diffuse insolation on a horizontal surface II dh is the difference between the global and 
the beam insolations. 

Solar radiation data is given for the locations of Viking Landers VL1 and VL2 based on the 
measured optical depths by the Viking Lander cameras. Solar radiation data at other latitudes are 
discussed in [3] and will be also presented elsewhere. The daily global insolation in k\V-hr/m 2 -day 
on a horizontal surface on Mars at VL1 and VL2 is shown in Fig. 4. The daily global insolation at 
VL1 and VL2 is repeated again in Figs. 5 (a) and (b) but together with the percentage of the diffuse 
insolation Hdh/Hh • The figure shows that during spring and summer where the optical depth is 
relatively low, the diffuse insolation comprises about 50 percent of the global insolation, whereas for 
autumn and winter, the optical depth is high and therefore the diffuse insolation comprises about 
90 percent. 

As mentioned before, there may be years on Mars which have no global storms (clear skies). 
The optical depth of the atmosphere for clear skies is typically 0.5. It is interesting to compare 
the insolation on the surface of Mars for clear skies with that for the opacities observed in 1977 by 
Viking [2,3]. This comparison is shown in Fig. 6. As one can see, the effect of the global storms is 
more pronounced at the lower latitude site VL1 than at the higher latitude site VL2. 


p (— t/ cos z)du> 


( 8 ) 


* These equations are in terms of actual (terrestrial) hours. A Martian day is 24.65 terrestrial hours long. In 
some cases it is convenient to define a “Mars hour” to be 1 / 24th of the Martian day, or 2.7% longer than a terrestrial 
hour. To get the insolation on Mars with reference to these Mars hours, replace the 24.65 in the numerator by 24 hr. 


4 


3.2.2 Inclined Surface 


The solar radiation on an inclined surface is based on solar radiation data on an horizontal surface 
Gh and Gdh* The global irradiance Gp on an inclined surface with an angle (3 is given by: 

Gp = Gb cos 9 + Gdh cos 2 (/?/2) + Gh sin 2 (/?/2) (9) 

where Gb is given in eq. (3) and 8 is the angle of incidence, the angle between the beam direction and 
the normal to the surface. The insolation on the surface is obtained by integrating the irradiance 
over the period from sunrise to sunset. The integration may be replaced by a summation, obtaining: 

Hp = £ Gbi cos 0,Aw + Y}Gdhi cos 2 (/?/2) + alGhi sin 2 (/?/2)]Aw (10) 

U/ C r W jr 

where the subscript i corresponds to the irradiance at solar angle and A oj is the solar angle 
interval. 

For collecting direct irradiance, the daily optimum tilt angle (3 for a south-facing surface in the 
northern hemisphere is: 

P = <f>-6 ( 11 ) 

where $ is the local latitude and <5 the angle of the sun from the Mars equator, which varies from 
+25° to -25° with the season, and /?= 0 corresponds to a horizontal surface. Averaged over a year, 
the theoretical optimal inclination angle for an array surface facing south is for the tilt to be equal 
to the latitude: ■ 


(3 = 4> ■ (12) 

The variation of the daily global insolati on Hp for a surface with j3 = <j> at VL1 and VL2 is 
shown in Fig. 7. The difference in the daily insolation for a horizontal surface and for the inclined 
surface with (3 — $ is shown in Fig. 8 as a percentage loss (or gain) compared to the horizontal 
surface. It is interesting to note that this tilt angle results in a decrease in the total yearly insolation, 
rather than a gain. One can observe that there is a loss in the insolation during the summer, when 
the optimum daily tilt angle (equation 11) is closer to horizontal than the (3 = 4> yearly optimum. 
While a surface tilted at the latitude angle is optimum for collection of the beam radiation, it 
is less efficient than a horizontal surface in collecting the indirect radiation. This decreases the 
performance of a tilted surface in winter and fall, when indirect radiation is the largest component 
of the insolation received. As a result, the advantage of a tilted surface is decreased considerably 
in all seasons by enough that it does not outperform the horizontal surface. The yearly loss of the 
insolation is small: about 2 percent at VL1 and about 1.5 percent at VL2. 

The insolation on an inclined surface /? = <f) for a Martian year with no global storms may 
be of interest. This is shown in Fig. 9 at VL1 and VL2 for a constant optical depth of r — 0.5 
throughout the year. 

3.2.3 Sun-tracking Surface 

Lastly, we calculate the insolation for a two-axis sun-tracking surface and compare the gain in 
insolation to horizontal and inclined surfaces. The global irradiance on a two-axis tracking surface 
is given by: 
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C*2x — Gb + Gdk cos 2 ( z / 2 ) + al Gh sin 2 (z/2) (13) 

and the daily insolation is obtained by integration or summation of the irradiance over the period 
from sunrise to sunset. 

The insolation on a two- axis tracking surface at VL1 and VL2 for a Martian year is shown in 
Fig. 10 for the observed opacities at VL1 and VL2, and in Fig. 11 for clear skies with r = 0.5. 
We also compared in Fig. 12 the daily global insolation on (1) a horizontal surface, (2) a fixed 
inclination surface with (3 = <j> and (3) a two-axis tracking surface for the observed opacities at VL1 
and VL2. 

The two-axis tracking surface results in a gain in total insolation received at both sites. The 
gain is +7 percent (compared to a horizontal surface) at the VL1 site, and +20.7 percent at the 
VL2 site. 

During the northern hemisphere spring and summer the insolation is relatively high, 2-5 kW- 
hr/m 2 -day, due to the low optical depth of the Martian atmosphere. These seasons, totalling a full 
terrestrial year, are the likely ones during which manned missions will be carried out. During the 
northern hemisphere autumn and winter, when dust storms typically take place, the insolation is 
considerably lower. 

4 Mars Environmental Considerations for Power Systems 

Other environmental factors of importance to photovoltaic system operation on the Mars surface 
are the temperature and the wind. The Martian surface temperature varies from a minimum of 
130° K, to as high as 300° K, with a mean of 215° K [6,7]. Air temperatures were measured by 
Viking at a height of 1.6 meters above the surface over a (Martian) year of measurement, including 
both local and global dust storms. Peak daytime temperatures varied from about 170° K at the 
VL2 site during a global storm, to almost 250° K near the summer solstice. Photovoltaic cell 
performance increases with decreasing temperatures, with peak efficiency occuring at 150-200° Iv; 
at lower temperatures the efficiency decreases. The temperature coefficient of efficiency depends 
on the material, and in general increases as the material bandgap decreases. Thus, low bandgap 
materials such as silicon and CuInSe 2 , which have high coefficients, increase in performance rapidly 
at the low temperatures to be found at Mars. 

Wind was also measured by Viking. The average wind speed measured at the VL2 site was 
about 2 m/sec [6], with winds of over 17 m/sec observed less than 1% of the time. On the average, 
the winds on Mars are not very high, and the extremely thin atmosphere of Mars means that the 
wind force is low. Abrasion by wind-blown dust of the same approximate composition as found on 
Mars has been measured on typical solar cell coverglass. Only wind velocities over 89 m/sec, higher 
than that measured by Viking, left abrasion on the coverglasses [7]. This finding is corroborated 
by the relative lack of abrasion damage on the Viking lander optical elements. 

A final aspect of the Mars environment not well characterized is the soil. One interpretation 
of the Viking lander life science experiments is that the Martian soil contains large amounts of 
peroxides and superoxides [6]. This hypothesis needs to be confirmed by direct chemical analysis 
of the soil. If it is confirmed, it will be important that solar arrays be built using materials not 
subject to attack by oxidants. 

There are several approaches to Mars photovoltaic power systems [9]. The first use will be to 
power unmanned rovers and scouting missions needed to investigate the geology (or areology) of 
Mars and scout out a location for the manned Mars base. For manned missions, the most important 
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criterion for a power system will be high reliability and absence of dangerous failure modes. Solar 
cell types to be considered include silicon, GaAs, and thin-film cells. Due to the high price of 
transporting materials to Mars, another priority for a surface power system is low weight. One way 
of reducing the weight of a photovoltaic array for Mars is to use extremely high efficiency solar cells, 
such as gallium-arsenide based tandem cells currently under development, with efficiency of 30% or 
higher. An alternative approach is to use thin-film solar cells, also now being researched. Thin-film 
cells have much lower efficiencies, but could be deposited in layers only a micron or two thick onto 
a lightweight plastic film which could be rolled out on the ground. Concentration systems, which 
can only take useful advantage of the direct component of the solar isolation, are less desirable on 
Mars, where direct insolation is only a fraction of the total energy. 

Like Earth, the duration of the night on Mars varies from winter to summer depending on 
the latitude. This is discussed further in [9]. An energy storage system is required for operation 
during the night, and in some cases storage may be desirable for continuous operation during local 
dust storms. An example power system including a storage system for night power is discussed 
by McKissock, Kohout and Schmitz [10]. Their power cycle was to provide 40 kW continuous 
power during the day, and a reduced power level of 20 kW during the night. For the night storage 
capability, they assumed that the hydrogen/oxygen regenerative fuel cell (RFC) has been developed 
to technology readiness. The RFC assumed pressurized gas storage. Cryogenic storage of reactants 
was determined to require too much equipment overhead to justify the slight advantages in storage 
volume. Round-trip efficiency for the storage system was 61.5%. 

Some resource processing technologies require thermal processing. One method of achieving high 
temperatures is the “solar furnace,” using concentrating mirrors or lenses. Since a solar furnace,' 
like a concentrator PV system, only utilizes the direct beam component of the solar insolation, the 
low direct beam insolation on Mars makes th e u tility of this approach much lower than in other 
environments. 

5 Conclusions 

The paper discusses solar radiation data on the surface of Mars from measurements made of the 
optical depths of the atmosphere at Viking Lander locations VL1 and VL2. 

The two global dust storms that were observed on Mars in 1977 makes it one of the worst years 
of dust activities. Even so, the opacities are relatively low during the northern spring and summer 
for more than a half Martian year (or more than a terrestrial year) for which the insolation is 
relatively high, varying from about 2 to 4 k\V-hr/m 2 -day on a horizontal surface, and varying from 
about 2.5 to 5.5 kW-hr/m 2 -day on a two-axis tracking surface. These seasons are the likely ones 
during which manned missions may be carried out. 

The Martian atmosphere contains suspended dust particles, resulting in significant scattering 
of light. This yields a large diffuse component of the solar insolation. This is presented in Fig. 
5. The diffuse insolation comprises about 50 percent of the global insolation during spring and 
summer and about 90 percent during autumn and winter. The effect of global dust storms on 
the insolation is shown in Fig. 6 where a substantial decrease in insolation is seen at VL1. The 
expected tilt at the yearly optimal angle /? = $ reduces, rather than increases, the insolation on 
the surface during spring and summer. The insolation during autumn and winter is affected by the 
high opacity rather than by the inclination of the surface. For a two-axis tracking surface, the gain 
in insolation is about 7 percent at VL1, and about 20 percent at VL2, compared to a horizontal 
surface. As mentioned before, the optical depth during the global dust storms dictates the amount 
of insolation on a surface independently of its inclination. 
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TABLE I - NORMALIZED NET FLUX FUNCTION COEFFICIENTS p(i, j, k) 


k -0 

H 

0 

1 

2 

3 

4 

5 

0 

1.002800 

-0.228681 

0.019613 

0.000231 

-0.000 1 30 

0.000003 

] 

-.450073 

1.335955 

-1.131691 

.402126 

-.063967 

.003758 

2 

5.566705 

-16.912405 

13.739701 

-4.756079 

,743740 

-.043159 

3 

-22.471579 

64.909973 

-52.509470 

17.997548 

-2 786548 

.160340 

4 

36.334497 

-101.800319 

79.895539 

-26.762885 

4.074117 

-.231476 

5 

-20.420490 

53 207148 

-39.949537 

12.977108 

-1.931 169 

.107837 

k - 1 

x 

0 

1 

2 

3 

4 

5 

0 

0.009814 

0 226139 

-0.1 17733 

0.030579 

-0.004090 

0.000218 

1 

-.156701 

.39682 1 

-.313648 

.099227 

-.013508 

.000651 

2 

1.161 122 

-3.758 1! 1 

3.007907 

-.987457 

.141693 

-.007320 

3 

-4 365924 

12.539251 

-10,394165 

3486452 

-.513123 

,027401 

4 

5.991693 

-17.498138 

14.191370 

4.765323 

.703675 

-.037960 

5 

-2.915099 

8.275686 

-6.593125 

2.173999 

-.320308 

.017335 
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Fig. 1: Optical depth at Viking Lander VL1. 
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Fig. 2: Optical depth at Viking Lander VL2. 
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4: Variation of daily global insolation II h ( kWhr/m 2 — day ) on a horizontal surface 
on Mars at VL1 and VL2 for a Martian year. 
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Fig. 6: Variation of daily global insolation Hh ( kWhr/rn 2 — day ) 
at VLl and VL2 for the observed opacities at VLl and VL2 
and for r = 0.5. 
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Fig. 7: Variation of daily global insolation Hp (kW hr / m? — day) on a tilted surface, /? = </>, 
at VL1 and VL2. 
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Fig. 8: Percentage of daily global insolation loss 
pared to a horizontal surface at VL1 and 
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Fig. 9: Variation of daily global insolation Hp ( kWhr/m 2 - day) on a tilted surface, (3 = <p, 
at VL1 and VL2 for r = 0.5. 
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Fig. 12: Variation of daily global insolation H (kWhr/m 7 - day) on a horizontal surface Hh , 
inclined surface Hp,(3 = <£; and on a two axis tracking surface for the observed 
opacities at VLl and VL2. 
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